R_n := 2 (Z^{3/2}/n^2) \sqrt{(n-l-1)!/(n+l)!} \exp(-Z r/n) (2Z/n)^l L^{2l+1}_{n-l-1}(2Z/n r).
The normalization is chosen such that the wavefunction \psi is given by \psi(n,l,r) = R_n Y_{lm}.